Conditions for Efimov Physics for Finite Range Potentials 
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We consider a system of three identical bosons near a Feshbach resonance in the universal regime 
with large scattering length usually described by model independent zero-range potentials. We 
employ the adiabatic hyperspherical approximation and derive the rigorous large-distance equation 
for the adiabatic potential for finite-range interactions. The effective range correction to the zero- 
range approximation must be supplemented by a new term of the same order. The non-adiabatic 
term can be decisive. Efimov physics is always confined to the range between effective range and 
scattering length. The analytical results agree with numerical calculations for realistic potentials. 

PACS numbers: 21.45.-v, 31.15.xj, 67.85.-d 



Introduction. Universal scaling properties in three- 
body systems arise when the scattering length a is much 
larger than the range ro of the underlying two-body po- 
tential [l[ . In this regime certain three-body observables 
are universal in the sense that they are model indepen- 
dent. This is colloquially referred to as Efimov physics 
@; 0, S, HI • Examples can be found in nuclear systems, 
small molecules, and particularly in cold atoms where the 
scattering length can be tuned to desired values using the 
Feshbach resonance technique. 

The universal scaling of Efimov trimers is usually said 
to exist for rms-sizes between r and a P, H, 0, @| • The ef- 
fective range R e from a low-energy phase shift expansion 
is sometimes used instead of r in this statement [a S Hi • 
This ambiguity occurs because ro and R e are often of the 
same order. However, for narrow Feshbach resonances in 
atomic gases R e can be much larger than ro @, and the 
implications for such systems need to be explored. 

Zero-range models, in particular in combination with 
the hyperspherical approximation 0, 0i have been 
successful in semi-quantitative descriptions of three-body 
systems in the universal regime. Semi-rigorous finite- 
range corrections have been attempted by includingthe 
higher order terms in the effective range expansion [a, 0] 
as a step towards the full finite-range calculations as in 



1 11 ] while maintaining the conceptual and technical 



simplicity of the zero-range approximation. 

The obvious generalization of the zero-range model is 
to substitute —1/a with — l/a+ (R e /2)k 2 , where k is the 
two-body wave-number, in the relevant expressions for 
the logarithmic derivative of the total wave-function at 
small separation of the particles. However, in three-body 
systems neither the two-body wave-number nor the small 
separation are uniquely defined, and rigorous inclusion 
of all terms of the given order is non-trivial. The lack of 
rigor in previous works could have serious implications 
for applications where finite-range effects are important, 
such as the stability conditions for condensates in traps, 
properties of cold atoms in lattices, and generally for Efi- 
mov physics. Experimental progress Q will soon require 



this increased accuracy near the boundaries of the uni- 
versal regime. 

In this Letter we derive, within the adiabatic hyper- 
spherical approximation [l2|, the rigorous asymptotic 
equation for the adiabatic potential, which includes the 
finite-range correction terms. The equation is suitable for 
the analytic studies of the finite-range corrections in the 
three-boson problem. We investigate the finite range cor- 
rections to the adiabatic potential and the non-adiabatic 
term and compare with the zero-range approximation. 

Adiabatic eigenvalue equation. We consider three 
identical bosons of mass m and coordinates interact- 
ing via a finite-range two-body potential V , where we 



assume V(rjk) = for r^u 



= TV 



- Tfc| > r . Only 



relative s-waves are included. We use the hyperradius 
r| 3 )2/x/3 and hyperangles tanevi = 



ik = \ri - (rj + r k )/2\ and \x 
[12]. In the following we shall 



P = ( r h + r h 
(i"jk/ri,jk)V3/2, where r ii3fi 
is an arbitrary parameter 
use one set of coordinates and omit the index. 

The adiabatic hyperspherical approximation treats the 
hyperradius p as a slow adiabatic variable and the hy- 
perangle a as the fast variable. The eigenvalue X(p) = 
v 2 {p) — 4 of the fast hyperangular motion for a fixed p 
serves as the adiabatic potential for the slow hyperradial 
motion. The eigenvalue is found by solving the Faddeev 
equation for fixed p > p c = 2ro^/J2, 



d 2 

da 2 



v 2 + U 



-2UK[<ip]. 



(1) 



Here tp(p, a) is the Faddeev hyperangular component, 

U{p,a) = V(psma/^Jl)mp 2 /(n 2 fi) (2) 



is the rescaled potential, and 
K[il>](p,a)=-= \ 

V3 J|f-a| 



ip{p,a')da' (3) 



is the operator that rotates a Faddeev component into 
another Jacobi system and projects it onto s-waves. The 



2 



total wave-function of the three-body system is ^(p, a) 
f(p)p~ 5 ^ 2 ®(p, a ) where 



$(p, a) 



V>(p,a) + 2^](p, a) 
sin 2a 



(4) 



The hyperradial function /(p) satisfies the ordinary hy- 
perradial equation [l^] with the effective potential 



V eS (p) 



m \ p z 



<*l|?l*>. (5) 



where Q is the non-adiabatic term and $ is normalized 
to unity for fixed p. 

We first divide the a-interval [0; tt/2] into two regions: 
(I) where U ^ 0, and (II) where U = 0. The regions are 
separated at a = ao where sinao = y/JIra/p = p c /(2p). 
In region (II) we have the free solution to Eq. (TTJ, 



( a ) = N(p) sm(va - v-) 



(G) 



with the boundary condition, ip 11 ^) = and normal- 
ization N(p). In region (I), since cto < 7r/6, Eq. (TTj) 
simplifies to 



da 2 



v 2 +U 



V/ = -2un[i! n ] : 



(7) 



with the solution ip 1 — t/j Ih — 2TZ[^p 11 ], where ip Ih and 
— 2TZ[4> 11 ] are homogeneous and inhomogeneous solu- 
tions, respectively. i[> Ih is the regular solution to 



m dr 2 



h 2 k 2 p 



V P (r) 



4> 



it, 



0. 



\fp 



(8) 
(9) 



Vm ' p 

where k p = ^fjlv / p and r — ap/^/Ji. When a — ► a^, 

oc sin(fc p r + S p ), (10) 

where the modified phase shift S p (k p ) arises from the 
modified two-body potential, V p . The solutions $ in re- 
gion (I) and (II) are now matched smoothly, leading to 



Aln^l =^-\nU H + 2TZ^ H })\ . 
da la ° da v L IJ]a o 



(11) 

After inserting Eqs. ([6]) and ([TO]) , this equation becomes 

8 



^ -!/cos(i>§) + ^sin(^f) 



Effective range expansion. In the limit p 3> p c , the 
p-dependent potential, V p (r), approaches V(r), and con- 
sequently S p approaches 5. The p-dependent low-energy 
effective range expansion corresponding to V p is then to 
second order 



k p cot S p (k p 



1 



Re(p) ,2 



a(p) 



(13) 



where a(p) and R e (p) are functions of 1/p 2 that converge 
to a and R e for p — > oo. Up to 1/p 2 in Eq. (TT3|) we get 



1 



a(p) 



P 2 



i?e(p) « i?e 



(14) 



The model dependent expansion parameter i?y , or "scat- 
tering length correction" , is found to be 



6ft 2 



U'(r)r 3 u(r) 2 dr, (15) 



where u is the zero-energy two-body radial wavefunction, 
asymptotically equal to 1 — r/a. Eq. ([12"]) then becomes 

^ -i/cos(i/f) + ^sin(i/%) _ 1 | fl ei /^ 

a 2 p 2 p 2 



P 



sin(i/f) 



(16) 

This equation without the last two finite range terms has 
the well-known purely imaginary solution = 1.00624i, 
or A = —5.0125, for p/y/pZ <C |a|. This solution gives 
Vcfi oc —1/p 2 which is the basis of Efimov physics. The 
i? e -term was included in [!, 0] , but not the model depen- 
dent iiV-term. The latter term makes the finite-range 
corrections to the zero-range adiabatic eigenvalues ex- 
plicitly non-universal. The last two terms in Eq. (|16|) 
restrict the solution A to the region |i? | ^ p/ \fp "C \a\, 
where 



(17) 



D ^ e 2 r> 

^0 — ~^ U Q ~~ H V> 



as seen in Fig. [T] where the lowest solution to Eq. (fTf 
is shown for different parameter choices. Thus, naively 
one would think that the lower limit for Efimov physics 
is determined by the model dependent length |-Ro|. How- 
ever, we will show later that Q restores universality and 
recovers the model independent effective range, R e . 

First, to illustrate the necessity of both l/p 2 -terms in 
Eq. (Tl7j|) we consider a large negative effective range cor- 
responding to a narrow Feshbach resonance 0] . To model 
the large |i? e | we pick an attractive potential with barrier 



P 



sin(i/|) 



k p cotSp(k p ), (12) V(r) =D sech 2 f x -) + Bexp ( - 2( x - - 2) 2 V (18) 



which defines v as function of p. The right-hand-side 
deviates from the zero-range approximations H, [l(| by 
using the rigorously defined phase shifts 5 P for V p instead 
of the original phase shifts 8. 



where D = -138.27, B = 128.49 in units of ft 2 /(mr 2 ), 
and x — 4.6667. The potential is negligible outside the 
range r$. The low-energy parameters are a = 556.88, 
R e = -142.86, R v = 73.031, and R = -0.71 in units of 
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FIG. 1: (color online) Adiabatic eigenvalues A(p) from 
Eq. (I16|l as function of hyperradius p, for large scattering 
length a and negative R e . Different values of the model de- 
pendent length i?o are used, showing that the universal solu- 
tion Ao exists in the region \Rq \ <C p/ \fp <C \a\. Lengths are 
in units of |7? e |. 
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FIG. 2: (color online) Exact numerical adiabatic eigenvalues 
A(p) for a potential with barrier, Eq. (|18[) (solid red line), 
compared to solutions of the eigenvalue equation, Eq. {IB) . 
The zero-range model (dotted line) includes only a. Crosses 
include a, _R e -terms and circles include a, R e , Rv-terms. The 
inset shows details around Ao- The arrow indicates the effect 
of the correction Ao — 2Ro /R e . 



ro. In Fig. [2] we compare X(p) obtained by exact numer- 
ical solution of the Schrodinger equation [ll[ containing 
the interaction Eq. (JTSJ) with the solution of Eq. (1 16[) . In 
the zero-range model (including only l/o), the —p 2 diver- 
gence for large p is below the numerical solution. At small 
distances, A approaches Ao, above the numerical solution. 
Inclusion of the i? e -term, as in [5, 7], provides a better 
large-distance behavior (since the dimer binding energy 
is corrected), but overshoots dramatically for p/yfp < a 
by approaching A ~ —4. Including consistently both 
R e - and -fiV-terms leads to complete numerical agree- 
ment with the exact numerical solution except for very 
small p- values where higher order terms are needed in 
Eq. (H 



Non- adiabatic corrections. We shall show that the 
non-adiabatic term restores model independence and re- 
covers |i? e | as the limit for the region of Efimov physics. 
For simplicity we only consider the limit \a\ = oo and 
assume |i?o| "C |-R e |- We first consider p/^pZ<^ \R e \ (re- 
gion A in Fig.[T]). Expansion of Eq. (fTB)) to first order in 
(y — vq) gives a small constant correction 



VqRc ^ ^ Re ^ 



(19) 



This correction is marked by the arrow in Fig. [2] (it is out 
of the range of Fig. Q}. This gives 



h 2 p / v%-l/4~2R /R e 
m v p 2 



Q 



(20) 



To evaluate Q we note that a large negative effective 
range (for |et| = oo) implies that the two-body wave- 
function u is localized mainly inside the potential range. 
Then the angular three-body wavefunction $ can be ap- 
proximated by uj sin(2a). The result is Q = cj p 2 , where 
c ~ —5/4 as confirmed numerically. This term cancels 
the main l/p 2 -part in Eq. (f20|) and hence prohibits Efi- 
mov physics for p/y/p <C |i? e |- The intuitive reason is 
that the two-body wavefunction is essentially zero out- 
side the potential, despite the large scattering length, and 
hence three particles can not interact at large distances. 
When p/yfp, » | .R e | (region B in Fig. [J) we find 



VqCq *—[l + 0{ — ) I, 

V / 



CO 



simVof )/vq 



^cos(^of) -cos(z/ §) + i/ f sin(^ f )' 



(21) 



(22) 



or Co — —0.671. This gives 

v M = #± (±^!± + S^E {RA _ 2Rv) _ 

m V p z p 6 

(23) 

The 1/p 3 dependence of the correction to the Efimov 
potential 1/p 2 was expected @. The model indepen- 
dent term proportional to R e /p 3 was recently calculated 
in [5j. However, we also get a model dependent term 
Rv/p 3 which is of the same order. Q generally receives 
contributions both from distances inside and outside the 
finite-range potential. Zero-range models only have the 
external part of the wavefunction, which depends on p 
only though the eigenvalue v(p). The zero-range result 
for Q is then 



Qzr 



M (^) 2 = M 0C0 V 2 ^, 



dp' 



(24) 



where M = ($\d 2 $/dv 2 )\ u=Vo . This fourth order correc- 



tion can be neglected in Eq. (|23p , as was done in [5j . How- 
ever, the internal part of the wavefunction contributes to 



4 



order 1/p 3 . To estimate this l/p 3 -term we take the ana- 
lytically solvable finite square well potential of range ro 
and \a\ — oo. This fixes R e = r and Ry = n 2 7r 2 ro/24 
where n is the number of bound states (including the 
zero-energy state) . We find 

Qb«« = <W(%-2iZv)^, ( 25 ) 
i p 3 

neglecting l/p 4 -terms. The model dependent Ry -terms 
in Eqs. (|23|) and (|25| cancel exactly, giving 

V^(p) = ^(^ + c ^-l)^\ (26) 
m V p z 2 p } / 

So the effective potential receives a R e /p 3 correction 
where the model dependent coefficient is different from 
zero-range models [5[ because of the inclusion of Q. We 
also expect the i?y-terms to cancel for general potentials. 
In conclusion, the Efimov effect persists for p/ ^fp \R e \. 

Atom-dimer potential. We have seen that model de- 
pendent corrections to Ao are cancelled by equivalent 
terms in Q. A similar effect occurs for the atom-dimer 
channel potential. Suppose the binding energy is Bo = 
H k D /m with corresponding wave number ko > 0. Then 
v = ikjjp/ JJL is an asymptotic solution to Eq. ()12|1 and 
A diverges as — p 2 corresponding to a bound dimer and 
a free particle. For this solution, the effective range 
expansion Eq. (fP5j) does not hold, since asymptotically 
k p — > iko is finite. Instead Eq. ([8]) reduces to the ra- 
dial two-body equation, with a normalized bound state 
s-wave function uo(r). Treating V„ — V oc 1/p 2 as a 
perturbation gives the correction 

A + 4 k 2 D f 0O 32 mV'(r) J 1 1 , , 0(T , 

—2- = -—-/ r u D—^r 1 &r—+0{ — ). 27 
P 2 P Jo P P 

Since lZ[ip] is exponentially small for the atom-dimer so- 
lution, Q can be computed using the unperturbed wave- 
function ip = y /pu D (ap/ y ^p), giving 

i r°° it 

Q = -T2+/ u D (ru' D +r 2 u' D )dr- + 0(-). (28) 
4p 2 J p 2 p 4 

By using the two-body radial equation and partial inte- 
gration the two integrals in Eqs. (|27f and (|28[) cancel. 
Thus the l//9 2 -terms in the effective potential Eq. ([5|) 
cancel exactly, giving V e g(p) = —Bp up to order 1/p 4 . 
Thus V c s only depends on R e through Bo- 
Effective range for Feshbach resonances. The effec- 
tive range near a Feshbach resonance has been estimated 
using a coupled-channels zero-range model Q , as 

R e = -2(AS?7iA/ia bs )" 1 , (29) 

where AB is the magnetic field width, Ap is the mag- 
netic moment difference between the channels, and a^g is 
the background scattering length. As an example we take 



the alkali atoms 39 K with the very narrow Feshbach res- 
onance at B = 825G having parameters AB = — 32mG, 
Ap = —3.92p B , and a^g = — 36a (HI- This gives the 
large effective range R e = —2.93 x 10 4 ao- For 39 K, r is of 
the order of the van der Waals length Z v dw = 1-29 x 10 2 ao 
4]. Since |i? e | 3> r , \R e \ determines the lower limit for 
Efimov physics and corrections to the universal regime 
are of order R e /a (not Z v dw/ a )- Thus, the window for 
universal physics is reduced. 

Summary and Conclusions. We consider a three- 
body system of identical bosons with large scattering 
length modelling a Feshbach resonance. The Efimov 
physics occurring in this "universal regime" is custom- 
arily accounted for by zero-range models. We use the 
adiabatic hyperspherical approximation and derive rigor- 
ously a transcendental equation to determine the asymp- 
totic adiabatic potential for a general finite-range poten- 
tial. We solve this equation for large scattering length, 
investigate finite range effects, and compare with exact 
numerical results. 

Inclusion of the effective range correction to the adi- 
abatic potential is insufficient in general. Crucial cor- 
rections of the same order must also be included from 
both the scattering length and the non-adiabatic term. 
These two contributions may separately be large but they 
tend to cancel each other. Accurate results in zero-range 
models must account for these new corrections. In con- 
clusion, the window for Efimov physics is precisely open 
between the effective range (not the potential range) and 
the scattering length. 
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